Two-dimensional, homogeneous, isotropic fluid turbulence with polymer additives 
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We carry out the most extensive and high-resolution direct numerical simulation, attempted 
so far, of homogeneous, isotropic turbulence in two-dimensional fluid films with air-drag-induced 
friction and with polymer additives. Our study reveals that the polymers (a) reduce the total fluid 
energy, enstrophy, and palinstrophy, (b) modify the fluid energy spectrum both in inverse- and 
forward-cascade regimes, (c) reduce small-scale intermittency, (d) suppress regions of large vorticity 
and strain rate, and (e) stretch in strain-dominated regions. We compare our results with earlier 
experimental studies and propose new experiments. 
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Polymer additives have remarkable effects on turbu- 
lent flows: in wall-bounded flows they lead to the dra- 
matic phenomenon of drag reduction [l|, ; in homoge- 
neous, isotropic turbulence they give rise to dissipation 
reduction, a modification of the energy spectrum, and 
a suppression of small-scale structures jsMll . These ef- 
fects have been studied principally in three-dimensional 
(3D) flows; their two-dimensional (2D) analogs have been 



studied only over the past decade in experiments 12. 13| 
on and direct numerical simulations (DNSs 
fluid films with polymer additives. It is important to 
investigate the differences, both qualitative and quanti- 
tative, between 2D and 3D fluid turbulence with poly- 
mer additives because the statistical properties of fluid 
turbulence in 2D and 3D are qualitatively different |17| : 
the inviscid, unforced 2D Navier-Stokes (NS) equation 
admits more conserved quantities than its 3D counter- 
part; one consequence of this is that the energy spec- 
trum for 2D fluid turbulence displays an inverse cascade 
of energy, from the energy-injection length scale linj to 
larger length scales, and a forward cascade of enstrophy, 
from linj to smaller length scales. We have, therefore, 
carried out the most extensive and high-resolution DNS 
study, attempted so far, of homogeneous, isotropic tur- 
bulence in the incompressible, 2D NS equation with air- 
drag-induccd friction and polymer additives, which we 
model by using the finitely-extensible-nonlinear-elastic- 
Peterlin (FENE-P) model for the polymer-conformation 
tensor. Our study yields several interesting results: We 
find that the inverse-cascade part of the energy spectrum 
in 2D fluid turbulence is suppressed by the addition of 
polymers. We show, for the first time, that the effect of 
polymers on the forward-cascade part of the fluid energy 
spectrum in 2D is similar to their effect on the energy 
spectrum in 3D fluid turbulence; in particular, there is a 
slight reduction at intermediate wave numbers and a sig- 
nificant enhancement of the fluid energy spectrum in the 
large- wave-number range; the high resolution of our sim- 
ulation is essential for resolving these spectral features. 



In addition, we find clear manifestations of dissipation- 
reduction-type phenomena j?!, S| : the addition of poly- 
mers to the turbulent 2D fluid leads to a reduction of 
the total energy, and energy- and enstrophy-dissipation 
rates, a suppression of small-scale intermittency, and a 
decrease in high-intensity vortical and strain-dominated 
regimes. We show that the probability distribution func- 
tions (PDFs) of cr^ and w^, the squares of the strain 
rate and the vorticity, respectively. a.gree qualitatively 
with those obtained in experiments [13| . We also present 
PDFs of the Okubo- Weiss parameter A = (w^ — cr^)/4, 
whose sign determines whether the flow in a given region 
is vortical or strain-dominated 18, 191, and PDFs of the 



polymer extension; and we show explicitly that polymers 
stretch preferentially in strain-dominated regions. 

The 2D incompressible NS and FENE-P equations can 
be written in terms of the stream-function ip and the 
vorticity a; = V x u(x, t), where u = {—dyip,dxtp) is the 
fluid velocity at the point x and time t, as follows: 
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Here Dt = dt + u.V, the uniform solvent density p — 1, 
a is the coefficient of friction, i/ the kinematic viscos- 
ity of the fluid, /i the viscosity parameter for the so- 
lute (FENE-P), and Tp the polymer relaxation time; to 
mimic soap-film experiments |13j we use a Kolmogorov- 
type forcing = kinjFo cos{kinjy)-, with amplitude 
-Fo; the energy- injection wave vector is kmj (the length 
scale linj = 27r/fci„j); the superscript T denotes a trans- 
pose, CjSj = {Rj3R-y) are the elements of the polymer- 
conformation tensor (angular brackets indicate an aver- 
age over polymer configurations), the identity tensor I 
has the elements Spj, fi^p) = (L^ — Tj/iL"^ — rp) is the 
FENE-P potential that ensures finite extensibility of the 
polymers, and rp = ■\/Tr(C) and L are, respectively, the 
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length and the maximal possible extension of the poly- 
mers; and c = fj./^v + fj,) is a dimcnsionless measure of 
the polymer concentration (20| . 

We have developed a parallel MPI code for our DNS, 
which uses periodic boundary conditions because we 
study homogeneous, isotropic, turbulence; our square 
simulation domain has side L = 27r and N'^ collo- 
cation points. We use a fourth-order, Runge-Kutta 
scheme, with time step St, for time marching and an 
explicit, fourth-order, central-finite-difFerence scheme in 
space and the Kurganov-Tadmor (KT) shock-capturing 
scheme [2l[ for the advection term in Eq. ([3]); the KT 
scheme (see Eq. (7) of Rcf. [§|) resolves sharp gradi- 
ents in Cp-y and thus minimizes dispersion errors, which 
increase with L and rp. We solve Eg. ^ in Fourier 
space by using the FFTW library [2^. We choose 
St ~ 5 X 10~* to 5 X 10~^ so that rp does not be- 
come larger than L. Table U lists the parameters we 
use. We preserve the symmetric-positive-definite (SPD) 
nature of C at all times by adaptingto 2D the Cholesky- 
decomposition scheme of Refs. H^l- In particular, 

we define J' = f{rp)C, so Eq. becomes 

DtJ = J.{Vu) + {Vuf.J - s{J qj, (4) 

where s = (L^ - 2 + f)/{TpL^), q = [d/{L^ - 2) - {L^ - 
2 + - 2)/(rpL2(L2 _ 2))], f Tt{J), and d = 

Tr[J.(Vu) + {Vvif.J]. Given that C and hence J are 
SPD matrices, we can write = CC^ , where £ is a 
lower-triangular matrix with elements tij , such that iij = 
for j > i; Eq.dH) now yields (1 < i < 2 and Pij = diUj) 
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EquationdS]) preserves the SPD nature of C if (a > 0, 
which we enforce 0, Q by considering the evolution of 
ln(£ii) instead of la. 

We use the following initial conditions (super- 
script 0): Cg^{x.) — Sp-y for all x; and w°(x) = 

^^p[sin(mi2; -|- ■m2y) + 



v[- A cos{Ay) + Ur^Y. 
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cos(mia; -|- TO22/)]"12/\/ ("t-i + m-i) ] • We maintain a con- 
stant energy- injection rate Einj =< Fu • u > by a rescal- 
ing of the amplitude Fq at every time step. Given a value 
of Einj , the system attains a nonequilibrium, statistically 



steady state after 
turnover time e 
square velocity. 



i 2Te — 3Te, whcrc the box-size eddy- 



In addition to a;(x, t), tp{-x.,t) and C(x, i) we 
obtain u(x, t), the fluid-energy spectrum E{k) = 
Efe-i/2<fc'<fc+i/2fc'^(IV'(k',t)P)t, where ()* denotes a 
time average over the statistically steady state, the total 
kinetic energy £{t) = (i|u(x, <)p)x, enstrophy Q{t) = 
(i|u;(x,OP>x, palinstrophy Vit) = (i| V x u;(x, t)^),, 
where ()x denotes an average over our simulation do- 
main, the PDF of scaled polymer extensions P{rp/L), 
the PDFs of w^, tr^, and = (w^ - cr^)/8, where 
= '^ij'^ij J 8'nd aij = diUj + djUi, the PDF of 
the Cartesian components of u, and the joint PDF of 
A and rp. We obtain the order-p, longitudinal, velocity 
structure function Sp{r) as follows: We first subtract the 
mean flow from the velocity field, u' = u — (u)t, then 
we define Sp{r) ^ {{[u' (r^ + r,t) - u'{r^,t)] ■ r/r}P)r„, 
where r has magnitude r and Tc is an origin, and 
denotes an average over time and the origin (we use 
= (*,j),2 < i,j < 5). We then extract the isotropic 
part Sp{r), by using an SO (2) decomposition [iMIiSl, by 
integration over the angle 9 that r makes with the x axis, 
i.e., Sp{r) = J^^ Sp{r)d9. We concentrate on 5*2(7') and 
the hyperflatness Fe{r) = SQ{r)/lS2{r)^]; the latter pro- 
vides a convenient measure of the intermittency at the 
scale r. 

In Fig. ([Iji) we show how the total kinetic energy 
£{t) (top panel), enstrophy il{t) (middle panel), and 
palinstrophy P{t) (bottom panel) of the fluid fluctuate 
about their mean values (£'(t))x, (ri(t))x, and (7-'(t))x as 
a function of time t for c = (pure fluid) and c = 0.4. 
Clearly (f (t))x, (r2(i))x, and (■P(t))x decrease as c in- 
creases. This suggests that the addition of polymers in- 
creases the effective viscosity of the solution. However, 
this naive conclusion has to be reflned because the ef- 
fective viscosity, arising from the polymers, depends on 
the length scale as can be seen by comparing the 

fluid-energy spectra, with and without polymers, in Fig. 
([Dd). Indeed, we can define the effective, scale- 

dependent viscosity Veik) = v + Ai^^k), with Av{k) = 
-MEA:-i/2<fe'<fe+i/2Uk' • iV-J).i.,/[Tpk^EP'"^ik)] and 
(V • J^)\i the Fourier transform of V • iJ. The right in- 
set of Fig. (Hi) shows that Au{k) > for fc < 30, 
where E^{k) < E-f{k), whereas, for large values of fc, 
Av{k) < 0, where EP{k) > E-^ {k); the superscripts / and 
p stand, respectively, for the fluid without and with poly- 
mers. The left inset of Fig. ([JJi) shows the suppression, 
by polymer additives, of n(fc) — T{k')dk' , where 

/ iZi(— k)Pij (k)(u X u))j(k)dil and the projector 
= Sij — -^^l^ , The suppression of the spectrum in 
the small-fc regime, which has also been seen in the exper- 
iments [lH and low-resolution DNS (see Fig. (4.12) of 
Ref. Q), signifies a reduction of the inverse cascade; the 
enhancement of the spectrum in the large-fc regime leads 
to the reduction in f2 and V on the addition of polymers 
(see the middle and bottom panels of Fig. (HJi)); to iden- 
tify this enhancement unambigouosly requires the high- 



r(fc) = 
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resolution run RIO, which is by far the highest-resolution 
DNS of Eqs. (U}-® (with 256 times more collocation 
points than, say, Ref. jl5l|). 

The second-order structure function S2{r) is related 
simply to the energy spectrum. It is natural to ask, there- 
fore, how S2{r) is modified by the addition of polymers. 
In Fig. (IH;), wc plot S2{r) versus r for c = (blue circles 
and run R7) and c = 0.2 (green asterisks and run R7); 
the dashed line, with slope 2, is shown to guide the eye; 
this slope is consistent with the <S'2('") ~ r'^ form that we 
expect, at small r, by Taylor expansion. At large values 
of r, S2{r) deviates from this behavior, more so for 
c = 0.2 than for c = 0; this is consistent with the results 
of the experiments of Ref. jT3|- The inset of Fig. ([1];), 
which contains plots of the hyperflatness ^6('') versus r 
for c = (blue circles) and c = 0.2 (green asterisks and 
run R8), shows that, on the addition of polymers, small- 
scale intermittency decreases as c increases. In Fig. dJK), 
we show how E'P(k) changes, as we increase c : at low 
and intermediate values of k (e.g., k = 1 and 30, re- 
spectively), EP{k) decreases as c increases; but, for large 
values of k (e.g., k = 100), it increases with c. Figure 
([2)d) shows how EP{k) changes, as we increase rp with 
c held fixed at 0.1. At low values of k (e.g., k = 1), 
E^{k) decreases as rp increases; but for large values of 
k (e.g., k = 100) it increases with rp. In Fig. ((^t) we 
give plots, for c = 0.1, of the spectra E'P{k) for L = 100 
(red triangles and run R8) and L = 10 (green asterisks 
and run R9); for comparison we also plot E^ [k) for c = 0; 
as L increases, the difference between EP{k) and E^ {k) 
increases at large values of k. 

We turn now to the PDF P{rp/L), which we plot ver- 
sus the scaled polymer extension rp/L in Fig. ([JJi) for 
c = 0.1 and L = 100 (red triangles and run R8), c — 0.4 
and L = 100 (black squares and run R8), and c = 0.1 * 
and L = 10 (green asterisks and run R9). These PDFs ^ 
fall off sharply (a) near rp/L ~ 1, because rp ^ L, and * 
(b) rp > V2, because rj, = Tr(C) > 2; at values of rp/L 
that lie between these two extremes, P{rp/L) shows a dis- 
tinct, power-law regime, with an exponent that depends 
on c, L, and Wi; as Wi increases, this power can go from 
a negative value to a positive value, thus signalling a coil- 
stretch transition. (We will give a detailed study of the 
L and rp dependence of P{rp/L) elsewhere.) Here we 
show a representative plot, for L = 6 and c = 0.1 (brown 
dots and run Rl), of P{rp/L) for a case in which the 
polymers are very close to a coil-stretch transition. 

In Figs. dSi), dSb), and ^) we show PDFs of A, 



and 



respectively, for c = (blue circles and 



run R7) and c = 0.2 (red triangles and run R7). These 
PDFs show that the addition of polymers suppresses 
large values of A, ct^, and uP' . If we make scaled plots 
of PDFs such as P{A/ Arms), then they fall on top 
of each other for different values of c; this also holds 
for P{cr^/crrms) and P{uj'^ /ujI^^). The inset of Fig. 
(|3^) shows a filled contour plot of the joint PDF of 



Tp and A; this illustrates that rp is large where A 
is large and negative, i.e., polymer stretching occurs 
predominantly in strain-dominated regions; this is 
evident very strikingly in Fig. ((SJi), which contains a 
superimposition of contours of rp on a pseudocolor plot 
of A; we give a video of a sequence of such plots at 
|http : / / www ■ yo utube . com/watch?v=5VGzQHiS8co&f eature=youtu 
and at http : //www. physics . iisc . ernet . in/~rahul/movies/polj 
The inset of Fig. (3c) shows that the PDF of any Carte- 
sian component of u is very close to a Gaussian. 

Our extensive, high-resolution DNS study of 2D fluid 
turbulence with polymer additives yields good qualita- 
tive agreement, in the low-fc regime, with the fluid-energy 
spectra of Ref. [l2[ , and the second-order velocity struc- 
ture function of Ref. 



13| . In addition, our study obtains 



new results and insights that will, we hope, stimulate 
new experiments. In particular, experiments should be 
able to measure (a) the reduction of (£'(<))x, (^^(0)x! and 
{'P{t))x shown in Figs. ([T^) and (IUd), (b) the modifica- 
tion of EP{k) at large k (Figs. ^) and (P)), (c) the c 
and L dependences oi E'P{k) illustrated in Figs.([2|a,)-([2l;), 
(d) the PDFs of {rp/L), A, cr^ and w^ (e) the stretch- 
ing of polymers in strain-dominated regions illustrated 
in Fig. dSJi), and (f) the suppression of the hyperflatness 
Ffi{r) at small r Fig. (HJ:). 
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St X lO" 
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V X 10" 


* Wi 


c 
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0.008 


10.0 


4.71 
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107, 85 


3.4, 3.6 
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100 1,2,4 


1.0 


0.005 
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2.26 4.52 9.04 
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221, 121, 53 , 38 
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1451, 1407 


8.0, 8.2 
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16384 
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0.2 


106, 61 
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TABLE I: The parameters for our DNS runs Rl-RlO for the two-dimensional, incompressible Navier-Stokes equation with air- 
drag-induced friction and polymer additives modelled via the FENE-P model. In all our runs the friction coefficient a = 0.01. 
We also carry out DNS studies of the two-dimensional NS equation with the same numerical resolutions as used in our runs 
with polymer additives. A'^'^ is the number of collocation points, St the time step, Einj the energy-injection rate, v the kinematic 
viscosity, and c the concentration parameter. The Taylor-microscale Reynolds number is Rex = ^TlEj yfvt and the Weissenberg 
number is Wi = rp \J d jv, where E is the total kinetic energy of the fluid and the energy dissipation rate per unit mass for 
the fluid. The dissipation scale is r\d = (v'^ /eY^'^ and kmax = A/3. 
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FIG. 1: (Color online) (a) Plots versus time t/T^ for run R7 of the total kinetic energy £ of the fluid (top panel), the enstrophy 
Q (middle panel), and the palinstrophy V (bottom panel) for c = (blue circles for run R7) and c = 0.4 (black squares for 
run R7); (b)Log-log (base 10) plots of the energy spectra E{k) versus k for c = 0.2 (red triangles for run RIO) and c — (blue 
circles for run RIO); left inset: energy flux n(A;) versus k for c = 0.2 (red line for run RIO) and c = (blue line for run RIO); 
right inset: polymer contribution to the scale-dependent viscosity Ai>{k) versus k for c — 0.2 (blue line for run RIO); (c) plots 
of the second-order velocity structure function 5*2 (f) versus r for c = (blue circle for run R7) and c = 0.2 (green asterisks for 
run R7); the dashed line with slope 2 is shown for comparison; the inset shows a plot of the hyperflatness Fe{r) versus for r for 
c = (blue circles for run R7) and c = 0.2 (green asterisks for run R7). 
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FIG. 2: (Color online) (a) Log-log (base 10) plots of the energy spectra E{k) versus k for c = (blue circles for run R6), for 
c = 0.2 (red triangles for run R6) and c = 0.4 (black squares for run R6); plots of E{k) versus c for Wi = 1.53 and k — 1 (left 
bottom inset), Wi = 1.53 and A: = 30 (left top inset), and Wi = 1.53 and k = 100 (right top inset); (b) log-log (base 10) plots 
from run R2 of E{k) versus k for Wi — 2.26 (blue circles for run R2) , for Wi = 4.52 (red triangles for run R2), and Wi — 9.04 
(black squares for run R2); plots of E{k) versus rp for c = 0.4 and k — 1 (left bottom inset) and c = 0.4 and k = 100 (right top 
inset); (c) log-log (base 10) plots, for c = 0.2 and Wi = 2.91, of E{k) versus k for L — 100 (red triangles for run R8) and L = 10 
(green asterisks for run R9); for comparison we show E{k) for c = (blue circles for run R8); inset: plots, for L = 100 and 
L = 10, of the scale-dependent correction to the viscosity Ai'{k) versus k (see text); (d) PDFs of the scaled polymer extensions 
P{rp/ L) versus rp /L for c = 0.1 and L = 100 (red triangles for run R8), c — 0.4 (black squares for run R8), c — 0.1 and L = 10 
(green asterisks for run R9), and c = 0.1 and L = 6 (brown dots for run Rl). 




FIG. 3: (Color online) Probability distribution function (PDFs) of (a) the Okubo-Weiss parameter A for run R7; the inset 
shows a filled contour plot of the joint PDF of A and rp, the square of the polymer extension, for run RIO, (b) a^, the square 
of the strain rate, for run R7, and (c) cj^, the square of the vorticity (the inset shows a PDF of the velocity component Ux for 
c = (blue circles for run R7) with a fit (1/2) exp(— u^/12.5) (blue solid line), and for c = 0.2 (green asterisks for run R7) with 
a fit (1/2.65) exp(— u^/20) (green solid line)); (d) a representative pseudocolor plot of A superimposed on a contour plot of rp 
for run RIO. 



